A unified model for vortex-string network evolution 
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We describe and numerically test the velocity-dependent one-scale (VOS) string evolution model, 
a simple analytic approach describing a string network with the averaged correlation length and 
velocity. We show that it accurately reproduces the large-scale behaviour (in particular the scaling 
laws) of numerical simulations of both Goto-Nambu and field theory string networks. We explicitly 
demonstrate the relation between the high-energy physics approach and the damped and non- 
relativistic limits which are relevant for condensed matter physics. We also reproduce experimental 
results in this context and show that the vortex-string density is significantly reduced by loop 
production, an effect not included in the usual 'coarse-grained' approach. 



INTRODUCTION 

Vortex-lines or topological strings can appear in a wide 
range of physical contexts, ranging from cosmic strings 
in the early universe to vortex-lines in superfluid helium 
(for reviews see ref. 0,0, El)- Gaining a quantitative un- 
derstanding of their important effects represents a sig- 
nificant challenge because of their nonlinear nature and 
interactions and because of the complexity of evolving 
string networks. Considerable reliance, therefore, has 
been placed on numerical simulations but unfortunately 
these turn out to be technically difficult and very com- 
putationally costly. This provides strong motivation for 
alternative analytic approaches, essentially abandoning 
the detailed 'statistical physics' of the string network to 
concentrate on its 'thermodynamics'. Here we present 
one such model for string network evolution, the velocity- 
dependent one-scale (VOS) model 0,0, and demonstrate 
its quantitative success by direct comparison with numer- 
ical simulations. We are able to describe the scaling laws 
and large-scale properties of string networks in both cos- 
mological and condensed matter settings. 

The first assumption in this analysis is to 'localise' 
the string so that we can treat it as a one-dimensional 
line-like object. This is clearly a good assumption for 
gauged strings, such as magnetic flux lines, but may seem 
more questionable for strings possessing long-range inter- 
actions, such as global strings or superfluid vortex lines. 
However, as we shall see, we will be able to establish good 
agreement between the VOS model and simulations in 
both 'local' and 'global' cases. The second step is to av- 
erage the microscopic string equations of motion to derive 
the key evolution equations for the average string veloc- 
ity v and correlation length L. This is a generalization 
of Kibble's original 'one-scale' model 0, and has been 
described elsewhere Iff. 

We make a detailed comparison between the VOS 
model and numerical simulations - currently the world's 
largest and highest resolution - using both direct field 



theory simulations of magnetic flux-lines, as well as sim- 
ulations of Nambu strings, treating them as localised line- 
like objects. We are able to demonstrate good agreement 
between all three approaches, thus underpinning the im- 
portant assumptions required for the VOS model. Fixing 
a single parameter, we are able to provide a good de- 
scription of cosmic strings throughout the history of the 
universe, from the friction-dominated regime after their 
formation, through the radiation-matter transition and 
into the accelerating epoch today. Significantly, we be- 
lieve the model also describes the same average features 
of an evolving vortex-line tangle in a condensed matter 
context, reproducing the expected and experimentally 
observed scaling law. The present VOS model in this 
context averages over both the background fluid friction 
and the Magnus force, but we believe it can be adapted 
further to incorporate other physical effects. 



THE VOS STRING NETWORK MODEL 

It is interesting to contrast the condensed matter and 
high-energy physics approaches to vortex string dynam- 
ics. In high-energy physics one uses, whenever possible, a 
one-dimensional view (obtained by integrating over radial 
modes of the vortex solution), described by the Nambu 
action. In this approximation the string equations of mo- 
tion are 0, 
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where x is the string position, e is a measure of the en- 
ergy along the string, and id = 2Tt + [x//3T 3 is a damping 
length scale; time derivatives are with respect to confor- 
mal time (for condensed matter purposes this is identical 
to physical time). All terms can be rigorously derived 
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from the Nambu action, except the last one in (JTJ which 
comes from the Kalb-Ramond action and describes the 
Magnus force arising when a global string moves through 
a relativistic superfluid background [8| whose density is 
parametrized by £ 0] • 

Now let us consider two limits of the equations of 
motion that are relevant in condensed matter systems. 
Firstly, if the damping term dominates, we find after 
some algebra (note that the Magnus force is not included, 
since it's not dissipative) 



(3) 



where the right-hand side is the friction force term, e.g., 
on a superfluid vortex H, 0, • Secondly, let us con- 
sider the non-relativistic limit but without damping. In 
this case we find 
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which is the equation describing the frictionless motion 
of a vortex filament in an unbounded fluid 9, 10]. 

By combining these two terms we can therefore repro- 
duce the equation of motion obtained by Schwarz |9j, llOl ] , 
which used an effective and more phenomenological ID 
approach (based on a coarse-grained order parameter) 
to obtain the terms one by one. Note however that the 
Schwarz equation has further (subdominant) terms, com- 
ing from non-local and boundary contributions (which we 
have neglected altogether). 

Now, it is well known experimentally 00 that vortices 
in such systems evolve following a scaling law L cx t 1 / 2 , 
and the same is true for example for defects in liquid 
crystals 0] . In fact this is expected to hold for any phase 
ordering with a non-conserved order parameter (and in 
this context it is sometimes referred to as the Lifshitz- 
Allen-Cahn equation). We will see below that this scaling 
law can be easily derived from an averaged version of 
the evolution equations (11 1211 — a process that might be 
compared to coarse-graining — in addition to also being 
relevant in high-energy contexts. 

The VOS model has been described in detail elsewhere 
0, 0], so here we limit ourselves to a brief summary. 
Any string network divides fairly neatly into two distinct 
populations, long (or 'infinite') strings and small closed 
loops. A phenomenological term must then be included 
to account for the loss of energy from long strings by the 
production of loops, which are much smaller than L — this 
is the 'loop chopping efficiency' parameter c. A further 
phenomenological term (characterized by a strength £ 
and a characteristic length scale Lj) is also included to 
account for radiation back-reaction effects. Note that all 
these parameters are constants — refer to 0, 0, ^| for de- 
tails. This term is negligible for GUT-scale Goto-Nambu 
string networks, but can be relevant in other cases. 



By suitably averaging Eqns. one can obtain the 

following evolution equations 

2^ = 2HL + cv + |V + 8E« 6 exp ( ~) , ( 5 ) 
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here H is the Hubble parameter (relevant for cosmology) 
and k is a velocity-dependent phenomenological (but oth- 
erwise universal) parameter, called the 'momentum pa- 
rameter', given by 0] 

^) = ^(l-, 2 )(l + 2V2,3)i^! ; (7) 

its non-relativistic limit is k nT (v) ~ 2\/2/ir. 

We can now take the condensed matter (non- 
relativistic) limit of the VOS model. All we need to do 
is set H = and the damping length to a constant. One 
finds a stable attractor solution 
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We have thus shown that our Goto-Nambu based mi- 
croscopic equations of motion and our averaged version of 
them successfully reproduce known condensed matter re- 
sults, respectively the Schwarz equation and the L cx i 1 / 2 
scaling law. Note that our solution demonstrates the im- 
portance of loop production — although this is not usu- 
ally included in theoretical or numerical analyses in the 
condensed matter context, it has been observed in exper- 
iments jl 1| . We will now further test our averaged model 
in the context of field theory simulations and cosmology. 



FIELD THEORY SIMULATIONS 

In a previous paper [l^ we have described the evo- 
lution of string networks in full 3D field theory simula- 
tions of the Abelian-Higgs model. This is a relativistic 
analogue of the Ginzburg-Landau theory of superconduc- 
tivity. Before reviewing the comparisons between these 
simulations and the VOS model, it is worthwhile noting 
a byproduct of this work which relates to vortex tangles 
in condensed matter. In order to create quiescent ini- 
tial conditions for string evolution in these simulations, 
instead of the relativistic equations, we began by solv- 
ing the corresponding diffusive equations (refer to [l^ 
for details). This was essentially equivalent to the non- 
relativistic evolution of magnetic flux-lines in a friction- 
dominated regime. Measurements of the string corre- 
lation length for the evolving network revealed a clear 
L oc t x l 2 scaling, as illustrated in Fig. ^ As well as con- 
firming the VOS model prediction in this case (JSJ, this 
has already been observed experimentally |ll| . 
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FIG. 1: Friction-dominated evolution of a vortex-string net- 
work. Measured values of the average correlation length dur- 
ing diffusive evolution of an Abelian-Higgs string simulation 
are plotted against the best-fit scaling law L oc t 1 ^ 2 
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FIG. 2: Fitting the VOS model to flat (top left), radiation 
(top right), matter (bottom left) and global (bottom right) 
field theory simulations, with c = 0.57, E = 0.5 (except that 
for global it's 1.1) and Ld = 47r. 



The relativistic evolution of the string networks also 
clearly revealed the predicted scaling laws and, remark- 
ably, the correlation length and velocities for all the sim- 
ulations had a good asymptotic fit from the VOS model 
using the single parameter c « 0.57. This fit was uni- 
versal regardless of whether the simulations were in flat 
space or in an expanding universe, or whether matter or 
radiation eras, as illustrated in fig. O There is no de- 
generacy here between the chopping efficiency c which 
determines the asymptotics and the massive radiation 
parameters S = 0.5 and Ld — 47r which only affect the 
initial conditions. However, for global strings with mass- 
less radiation Ld — > oo there is such a degeneracy between 
c and X because they act in the same manner asymptot- 
ically. However, assuming the same loop chopping ef- 
ficiency c = 0.57 for local and global cases, requires a 
much higher damping coefficient £ = 1.1 for the latter 
as expected for massless radiation [l2|. (These results 
and fits are also in reasonable agreement with other re- 
cent simulations of global strings in ref. (l3j-) This ex- 
cellent correspondence for both local and global strings 
appears to establish the validity of the two key ('localiza- 
tion' and 'thermodynamic') assumptions underlying the 
VOS model. 



HIGH-RESOLUTION GAUGE STRING 
SIMULATIONS 

We now compare our model with ultra-high resolu- 
tion numerical simulations, the details of which will be 
reported elsewhere. Specifically, we discuss both the 
approach to the linear scaling regime in the radiation 
and matter epochs, and the transition between the two 
epochs. Note that in this case both friction and radiative 



back-reaction are negligible. We first compare the model 
with three different ultra-high resolution runs (75 points 
per correlation length, fixed physical resolution) in the 
radiation and matter epochs, see Fig. [31 For all 6 runs, 
the analytic model fit displayed corresponds to a loop 
chopping efficiency parameter c = 0.23. We find that this 
value also approximately reproduces the earlier results of 
Bennett and Bouchet [yj and Allen and Shellard [lfij . 
Residual differences can be attributed to the much higher 
resolution of our runs. Earlier work [l6| reported differ- 
ent scaling values, but numerical codes were substantially 
improved subsequently. Thus fixing this parameter via 
the radiation era, our model correctly predicts the mat- 
ter era scaling large-scale properties without any further 
tampering with parameters. We estimate the loop chop- 
ping efficiency to have the value c ren = 0.23 ± 0.04. We 
emphasize that we expect this to be a 'universal' param- 
eter, independent of the cosmological scenario in which 
the string network is evolving. However, if one performs 
analogous simulations in flat (Minkowski) spacetime, one 
does find a different value, c ren = 0.57 ± 0.04, which co- 
incides with the value we found above for field theory 
simulations. This is because the amount of small scale 
structure present in Goto-Nambu expanding runs is much 
larger than that in field theory and/or Minkowski space 
runs, and this has an influence on the large-scale features 
of the network described by the model. Hence the two 
values can be regarded as the 'renormalized' and 'bare' 
chopping efficiencies. This is discussed at greater length 

Now, given that we can reproduce the radiation and 
matter era scaling values, can we also reproduce the tran- 
sition between them? Note that this point is not straight- 
forward: nothing guarantees us, a priori, that the model 
will produce the correct timescale for the transition even 
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FIG. 3: The approach to scaling of radiation (left) and matter 
(right) era ultra-high resolution simulations. Solid lines are 
data from three different runs (for each epoch). All dashed 
lines are solutions of the analytic model, with the relevant 
initial conditions, for c = 0.23. 

if both its asymptotic limits are correct. Here it is not 
possible to do a single numerical simulation spanning all 
the relevant time interval, so we have performed a se- 
quence of twelve different high-resolution runs (16 points 
per correlation length, fixed physical resolution) which 
together span the required range, see Fig. Again we 
can see that the analytic model, with the previously de- 
termined parameter c — 0.23 can actually do a good job 
in predicting the overall timescale of the transition be- 
tween the regimes. We also emphasise that a GUT-scale 
network at the present time has not quite yet reached 
the asymptotic matter-era scaling regime. This slow re- 
laxation process has been pointed out before and is 
clearly of vital importance for an adequate analysis of 
cosmic string structure formation scenarios. 

DISCUSSION AND CONCLUSIONS 

We have described a simple analytic model for vortex- 
string evolution and shown how it successfully reproduces 
the large-scale features of numerical simulations of both 
Goto-Nambu and field theory string networks, as well as 
of experiments in condensed matter physics. This quan- 
titative correspondence provides strong evidence in sup- 
port of the main assumptions on which the VOS model is 
based, notably string 'locality' and 'thermodynamic' av- 
eraging. Our results confirm that the dominant mecha- 
nisms affecting string network evolution are loop produc- 
tion and damping, whether from friction or radiation de- 
pending on the context. In condensed matter systems we 
find that the two are of comparable magnitude, despite 
loop production being neglected in the usual treatments. 
For global strings, loop production is dominant but radia- 
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FIG. 4: The evolution of a GUT-scale string network in the 
transition from the radiation to the matter-dominated epochs. 
Solid lines are data from a set of 12 different high-resolution 
runs. The dashed curve is the solution of the VOS model with 
c = 0.23, while the dotted curves correspond to c = 0.21 and 
c = 0.25. 



tive damping can significantly affect the network density. 
The outstanding issue raised by comparisons of gauged 
string networks in Nambu and field theory simulations 
remains the modelling of small-scale features. A more de- 
tailed analysis of the simulations provides some valuable 
insights into this problem, but this will be the subject of 
a forthcoming publication. 
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